Entanglement spectra of complex paired superfluids 
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We study the entanglement in various fully-gapped complex paired states of fermions in two 
dimensions, focusing on the entanglement spectrum (ES), and using the Bardeen-Cooper-Schrieffer 
(BCS) form of the ground state wavefunction on a cylinder. Certain forms of the pairing functions 
allow a simple and explicit exact solution for the ES. In the weak-pairing phase of ^-wave paired 
spinless fermions (I odd), the universal low-lying part of the ES consists of \l\ chiral Majorana 
fermion modes [or 2\i\ (I even) for spin-singlet states]. For \l\ > 1, the pseudo-energies of the 
modes are split in general, but for all I there is a zero-pseudo-energy mode at zero wavevector if 
the number of modes is odd. This ES agrees with the perturbed conformal field theory of the edge 
excitations. For more general BCS states, we show how the entanglement gap diverges as a model 
pairing function is approached. 



PACS numbers: 74.20.Fg, 74.20.Rp, 03.67.Mn 



The hunt for both theoretical and experimental meth- 
ods that would allow to distinguish between two dif- 
ferent topological phases of matter [l], [2} has been re- 
cently revived, greatly due to the question whether the 
experimentally-observed fractional quantum Hall effect 
(FQHE) at v = 5/2 supports excitations with non- 
Abelian statistics 0, HI ■ An experimental observation 
of these non-Abelian anyons, possibly in interferome- 
try experiments 043]) would be a major breakthrough 
and an important step forward in the field of topological 
quantum computation 8]. Other proposals for experi- 
mental realizations of non-Abelian phases of matter in- 
clude two-dimensional (2D) p-wave superfluids [9], such 
as films of some phase of 3 He [l(|, ultracold Fermi gases 



(atoms interacting via p-wave Feschbach resonance 11 1 
or micro- wave dressed polar molecules fl2j ) , or closely re- 
lated physics at the interface between a superconductor 
and a topological insulator [3] ■ 

On the theoretical side, our understanding of the col- 
lective behavior of electrons and their emergent prop- 
erties (such as properties of the excitations, including 
fractional char ge a nd statistics) relies mostly on trial 



wave- functions [14l Il5j . Numerical comparisons between 
the trial wave-functions and the realistic ones (e.g. for 
Coulomb interactions) are most meaningful when they 
focus on quantities that are robust inside a topological 
phase. In recent years, quantum information concepts 
[l6l [13] such as quantum entanglement have provided 
valuable insights in this matter. First, the topological en- 
tanglement entropy (EE) associated with a reduced den- 
sity matrix of a 2D ground state in a topological phase 
la . Il9j was shown to take a constant value throughout 
the phase. Two different phases, however, might have 
the same topological EE. More refined approaches con- 
sider the full structure of the reduced density matrix. Its 
set of eigenvalues was studied numerically by Li and Hal- 
dane (LH) [2c| . They defined the entanglement spectrum 
(ES) as the set of real pseudo-energies e n appearing in 



the Schmidt decomposition of the ground state for a 
(spatial) bipartition of the Hilbert space rl = Ha ®T~Lb'- 

1 



e £ " /2 \-4>A,n) ® |^fl,n) 



(1) 



with Z = E„ e_£ "- For the FQHE at v = 5/2, LH ar- 
gued that as the size goes to infinity, the ES contains 
a universal low-pseudo-energy part (i.e., separated by a 
gap from the rest of the ES), which is related to the physi- 
cal edge spectrum of this phase [IjJ (such a relation with 
a chiral conformal field theory was also anticipated in 
Ref. [Hj]). This entanglement gap varies inside the topo- 
logical phase; it goes to infinity for a particular model 
wave function: the Moore- Read (MR) wave function 

The ES has been studied in various systems, includ- 
ing FQHE 22, 23 1, spin chains or ladders [24| and lattice 
complex paired superfluids [25|. Most of the results rely 
on numerical calculations of the ES. Exact analytical re- 
sults in basic examples or toy models would shed light 
on this topic. However, such progress has been largely 
precluded by the intrinsic technical difficulty of studying 
these highly entangled phases of matter. 

In this Letter we study 2D complex paired superfluids 
starting from BCS theory 14], and show that the ES 
can be worked out explicitly for certain model wave func- 
tions. The problem essentially reduces to a free fermion 
problem, for which various other methods are known [26| . 
To the best of our knowledge, these tools have not been 
used so far to make analytic calculations of the ES in 
topological superfluids, however, there are general argu- 
ments about the relation with the edge spectrum [27H29I . 

Complex paired superfluids — In the BCS theory [14j 
for translation-invariant systems, the fermions in the 
ground state form pairs, the members of which carry op- 
posite momenta k and — k. The mean-field Hamiltonian 
for spinless or spin-polarized particles is 



Hbcs = 



^kcj^Ck + i ^A k c_ k c k + A k c k cL 



(2) 



2 



,2 

where £k = — M is the kinetic energy minus the chem- 
ical potential for each single particle, Ak is the gap func- 
tion (Ak its complex conjugate), and cj^Ck are fermion 
creation/annihilation modes. The ground-state of @ is 



where all the particles have transverse momentum ±k y 
and are within part A (x < 0), and a similar definition 
holds for |f2s) fe with x > 0. Then, the (normalized) 
wave function \ fl) k is created from the cut ground state 



|f2) oc exp 



2 

k 



|0), 



(3) 



5k = (£k — V^k" HAkPJ /^k is the -pairing function. 

Initially we will assume that, as a function of k, Ak is 
an eigenfunction of rotations, with angular momentum 
i (the relative angular momentum of the particles in a 
pair). The Fermi statistics of the particles allows only 
odd values of l\ in such states, parity and time-reversal 
symmetry are broken. We first consider I = — 1 (I = +1 
is similar). In two dimensions, the I = — 1 gap func- 
tion behaves generically as Ak — A(k x — iky) at small 
k (A > without loss of generality). At small k, the 
behavior of <?k depends on the sign of u, and there is a 
topological phase transition at /i — [9]. In the strong- 
pairing phase (/x < 0), ~ (k x — ik y ) as k — > 0, and 
the pairing function in real space decays exponentially 
with the distance between the particles in a pair. At 
large scale, the fluid can be thought of as a Bose conden- 
sate of point-like Cooper pairs. In the weak-pairing phase 
(p. > 0), gk ~ l/(fcx + iky), and the pairing function de- 
cays algebraically. At large distances, the component of 
the ground state ((3]) with N particles at complex coordi- 
nates Zj — Xj + iyj (j = 1, . . . , N) is given by @ 



tp(zx, . . .,z N ) 



|fi) - Pf 



(4) 



where Pf denotes the Pfaffian of a matrix. This Pfaffian 
form is a factor in the MR wave function [2| , and its long- 
range behavior is the first indication of the topologically 
non-trivial nature of the weak-pairing p ±ip phase Q . 

Schmidt decomposition on the cylinder — We consider 
the BCS wave function (j3]) on an infinite cylinder of cir- 
cumference L y , with coordinates (x,y) € Kx [0,L y ). Ro- 
tational invariance of the cylinder implies that the ground 
state is a product: 



(5) 



fc„>0 



where the ±k y 's are the discrete y-components of mo- 
mentum for the particles (k y L y /2TT is either integer or 
half-integer depending on the boundary conditions of the 
fcrmions). The Schmidt decomposition for a transverse 
cut of the cylinder at x = (Fig. [TJ is obtained as fol- 
lows. Let us call A (resp. B) the semi-infinite cylinder 
with x < (resp. x > 0). First, one introduces the 



(normalized) cut ground state \£Ia) 



k y > 0, with 



in, 



oc e 



/-oo dx + dx - 9k y (X-,X+)cl_ _ 



ns). for each 



«y\0), (6) 



\n) k cx |0 A ), ® |n B ) fc , 



(7) 



with a sharing operator G\ generating shared pairs (one 
particle in A, one in i?) 



/OO />OC 
dx + / dx_ 5 ^(a ; _,x + )cL,_ 
-OO J — OO 



where g ky B (x-,x + ) = g ky (x-,x + ) if (x_,x + ) eixB 

or (x_,a; + ) e B x A, and otherwise. (i.e. ) can 
be Schmidt decomposed 



(9) 



where fe n (d B k n ) creates a particle in one of an 
orthonormal set of single-particle states lying entirely in 
A (resp., £?), r is the Schmidt rank of G\ , and a n (k y ) 
are some coefficients. Expanding the exponential in l[7]). 
one thus gets 2 r terms. Generically, 2 r is the size of the 
smaller Hilbert space T-La or Hb, and is infinite in our 
continuum BCS model, so it is difficult to make progress 
from eq. ([7])- Moreover, the different terms in the expan- 
sion are not orthogonal in general. 

FIG. 1: (color online) Origin of the two terms in the Schmidt 
decomposition (|13[) . (a) Independent configurations of pairs 
in A and B. (b) Shared Cooper pair between A and B. 



Model BCS wave functions — However, for one trial 
form an exact solution is possible. This is the case when 
Ak = A(k x — iky) for all k, and fi — m*A 2 /2, then 



i* A 



3k 



(10) 



for all k, and in the plane g(x,y) oc 1/z for all x, y. (This 
special point has been studied in Ref. [30l|.) For a finite 
density of particles, an ultraviolet cutoff is needed, which 
for the cylinder we take to be a cutoff at large k y . The 
pairing function in terms of k y > and x± , eq. (| 10|) , is 



9k.. 



.(x-,x + ) = im*Ae- k « (x -- x+) Q(x- -x+), (11) 



where <d(x) is the Heaviside step function, and so the 
shared pairing function g k B (x- , x+) oc e~ kyX ~ e kyX+ is 



3 



factorized (x+ must be in A and X- in B). This leads to 
a Schmidt decomposition for Gj. of rank r = 1: 



G. = 



(12) 




Thus, expanding the exponential in (0, we get only two 
terms. The first term corresponds to configurations with 
independent Cooper pairs in A and B, and the second one 
to configurations with one Cooper pair shared between A 
and B (Fig. Q}. The two terms are obviously orthogonal, 
because (for k y ^ 0) they have different total transverse 
momentum in A (or B). Their norms can be computed 
using Wick's theorem and two-point correlation functions 

in the original ground state |fl): (d,Ad A ^ — (dsd^ = 
vj\J\ + v 2 and (d^d^ = v (l — u/yT + v 2 ) , where v — 
k y / (m* A). The Schmidt decomposition of |f2) fc becomes 



l*V = 



\/T+ e" £ ( fe a) 



in 



+ e - e (*„)/a b t 



with 



(13) 



\Qa) and 6 



B,fc„ 



k v \^a), normalized such that (£Ia\ bA,k y bA k 
(CIb\ bBfkybjj j. |S1b) = 1, and the pseudo-energy is 



|fU) oc 

\n A ) = 



e(fc v ) = 2 In \J1+ (k y /m*A) 2 + k y /m*A 



(14) 



Chiral ES— Thus, the ES of |fi) for the model state 
consists of a fermion mode of pseudo-energy e{k y ) for 
each ky > appearing in ([5]). If we introduce fermion 
operators rf k = r;_ fc for all k = ±k y , = h+k',o, 

then the spectrum of the pseudo-Hamiltonian 



H 



ES 



KM V-kyVky 



reproduces the ES. At small k y , e(k y ) ~ 2k y /(m*A), so 
the ES is that of a chiral Majorana fermion theory. In the 
case of periodic boundary condition in the y-direction, 
ky = is included in the product ([5]). The expressions 
in eq. (H"2"j) are divergent, and the x integrals should be 
cutoff at x = ±L x /2. There are again two terms in the 
Schmidt decomposition of |fi) fc =0 , and e(0) = 0. We 
have incorporated this above by introducing a self-adjoint 
("Majorana") operator 77J = ijo (which requires a two- 
dimensional space of states for k y = 0). Another such 
operator (located far away along the cylinder) would be 



needed to make a complex fermion operator Jsimilar to 
the case of the edge spectrum on a cylinder [211] V 

Higher angular momenta — The generalization to 
gaircd states with general £ < — modeled by Ay = 
A(k x —ik y )^ in the Hamiltonian (0 — reveals an interest- 
ing structure. For odd angular momenta £ in the weak- 
pairing phase /i > 0, we consider the model pairing func- 
tion = —2fj,/[A(k x + ik y )W] for spinless fermions. For 
even £ we consider the spin-singlet case with the same 
form of model pairing function (appearing in the wave 
function as <?kclk± c kt)- I* 1 contrast to the I = — 1 case, 
there is no value of \i for which these are the exact func- 
tions, but they still capture the correct long-distance be- 
havior of the weak-pairing phase. We find (for fixed spin 
orientation for even £) 

g ky (x-,x+) oc {x- - x+f^e-^-^Oix- - x+). 

The sharing operator ^ has rank r — \£\, leading to a 
Schmidt decomposition of |£!) fc of rank 2^. The cal- 
culation of the pseudo-energies can be done as in the 
p-wave case using Wick's theorem, although this now in- 
volves diagonalizing a matrix of rank \£\. For example, 
for £ = — 2, recombining the two spin sectors (t or 4- for 
particles with momentum +k y ), we find that the ES at 
small k y corresponds to 



H E s = 



k v >0 L 



y/2Ak y 



Va.s^ — ky'Ha^s.ky ~t~ M a bT] a s fc y 7/6,s ,k y 



with sums over repeated indices a, b = 1, 2 and s =j, 
4_. Here M — ma y is a splitting matrix 31], and 
m = 21n(l + v2). The pseudo-energy eigenvalues are 
then given by e±(k y ) = ±m + [y/2Ak v / fx) and are spin- 
degenerate. We expect a similar form in general, with 
\£\ chiral Majorana fields (2\£\ with spin) and an \£\ x \£\ 
antisymmetric Hermitian splitting matrix M (the term 
linear in k y might contain different velocities for differ- 
ent fermion types). Thus when \£\ is even, the modes at 
ky — > come in pairs with pseudo-energy ±m a , a = 1, 
. . . , \£\/2 (with spin degeneracy). When \£\ is odd, M has 
at least one zero eigenvalue, and generically we are again 





(15) left with one Majorana zero mode operator at k y 



only. This is consistent with arguments that these weak- 
pairing phases of £-wave superfluids should support non- 
Abelian vortex excitations when £ is odd, but not when 
t is even [9]. (Another case is spin-triplet odd-£ pairing, 
which for model pairing functions like those above leads 
to 2\£\ modes, with £ odd 0.) 

We remark that the model wave functions we have 
used possess an interesting property, that is directly re- 
lated to the number of Majorana fields which appear in 
the ES. The configurations with non-zero amplitude in 
the wave function of |S!) fe (k y > 0) obey some specific 
rules. For \£\ = 1, the fermions with +k y and with — k y 
must alternate in x (Fig. [IJ) . Similarly, for \£\ = 2, there 



4 



can be at most two successive fermions with identical 
momentum-spin +k y but not three, and the same for 
— k y 1 fermions. For general I, configurations with up to 
\£\ successive +k y (resp., —k y ) momenta can occur, but 
not more. Similarly to the p-wave case (Fig. [T]), there 
are several possibilities for the configurations of particles 
across the cut between A and B, leading to the 2' e ' (or 
2 2 ^') terms in the Schmidt decomposition (|f 3|) . 

Entanglement gap — In the weak-pairing phase, for 
7^ m*A 2 /2, the pairing function can be written 



,9k = /(|k| 2 ) 5k 0) with /(|k| 2 ) -> 1 when |k| -> 0, 



and 
sider £ 



is the model pairing function (|10j) 
-1 for definiteness). For fixed 

there is a new length scale ro(k y ) = 



con- 



(we 
k v > 0, 



the 



_ \f'(k 2 y)\ 

one-dimensional wave function |f2) fc . The pairing func- 
tion in real x-space g ky is a convolution of 
with a kernel of width ~ ro(k y ). In other words, for 
x+\ > r (k y ), one has g ky (x-,x + ) ~ g ( k °\x-,x + ), 



\x- -x+\ ^> r [K y ), one nas g ky (x-,x + ) ~ g ky '(X-,x A 
but the behavior of the pairing function is different for 
\x- — x + \ < ro(fcy). This produces (infinitely many) 
additional terms in the decomposition of the correspond- 
ing sharing operator G ky The additional terms in the 
Schmidt decomposition come from configurations with 
pairs of particles with momenta ±k y in the region close 
to the cut —r (k y ) < X- < 0, < x + < fo(k y ) 
(or with X- and x+ exchanged). Close to the special 
point [i — m*A 2 /2, ro(k y ) is small, and one expects 
such terms to be of order ~ r (k y ) 2 (p±k y (~r )p Tky (r )) 

where p ky (x) = c x k c Xj k y is the density of particles with 
momentum k y . Apart from a small change in the chiral 
mode (TT4")) (which must still behave linearly when k y — > 0, 
with a coefficient depending on /x), the only effect is that 
additional pseudo-energies appear in the ES, of order 
-In [r {k y ) 2 (p±k v {-r )p T k v {rQ))]- At k y = we get 



a rough estimate of the gap ~ — In 



1 - /i/(m* A 2 /2) 



Strong-pairing phase — In the strong-pairing phase, 
the exponential decay of the pairing function may be 
approximated by (derivatives of) a delta function, e.g. 
9k y (x) oc (d x + k y )S(x) for I = — 1. The term k y 5(x) is 
local, so it does not generate any entanglement. The term 
5'(x) might give some non-trivial contribution to the ES, 
but one can repeat similar arguments as for the entangle- 
ment gap. A regularization of 5'(x) would be localized 
on a short length r$. For a finite density of particles 
(p{k y )), one would then find that all the pseudo-energies 
go to infinity as ro — > 0. We conclude that, provided we 
keep the density of particles (p(k y )) finite, the ES in the 
strong-pairing phase cannot contain a universal gapless 
part. This is expected, because the strong-pairing phase 
is adiabatically connected to the vacuum as /i — > — oo 
@. Note that for s-wave pairing, I = 0, there is only the 
strong-pairing phase. 

Discussion — In more general states of BCS form, the 



rotational covariance of Ak may be broken, as can spin- 
rotation symmetry, and even translation symmetry. Be- 
cause of the chiral form of the low-lying part of the ES 
(and assuming a local Hamiltonian), the total number 
of chiral Majorana modes in the ES found here must be 
robust under such perturbations, unless the system un- 
dergoes a phase transition at which the bulk energy gap 
vanishes. Our approach can also be applied to paired 
states of other symmetry classes (including those with 
time- reversal symmetry), to multi-component fermions, 
and in higher dimensions (32[. 

In conclusion, we have shown that for certain model 
wave- functions for the complex ^-wave superfluids the ES 
can be computed exactly, and the entanglement gap is 
infinite. For spinless £-wave {I odd) pairing, the ES con- 
sists of |£| chiral Majorana fields (2|£| for spin singlet). 
The low-pseudo-energy part of the ES is that of a per- 
turbed chiral CFT which is also that of the edge states 
13, SI SHI, and is expected to be universal. 
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